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Abstract 

We discuss various contributions to the event-by-event charge-asymmetry dependence of 7r + 
and tt~ elliptic flow, recently measured by the STAR Collaboration at RHIC. It is shown that 
under general assumptions, the difference between and at a given fluctuating value of an 
asymmetry parameter, A, is a linear function of A, as observed in the preliminary data. We discuss 
two mechanisms that are qualitatively consistent with the experimental data and result in a signal 
of the correct order of magnitude. Our subsequent hydrodynamic calculations, assuming local 
charge conservation at freeze-out, yield a qualitative and partial quantitative understanding of the 
observed signal, so offering a detailed test of the hydrodynamic model in heavy-ion collisions. 
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I. INTRODUCTION 



Recent experimental measurements of elliptic anisotropy have challenged the conventional 
hydro dynamic picture of relativistic heavy-ion collisions that had proven very effective in 
facilitating our understanding of many features of particle production in such collisions [1-3]. 

For instance, as reported in Ref. [4], at a given energy and centrality of a collision the 
elliptic flow for particles and antiparticles is slightly different. In particular the elliptic flow 
for protons is larger than for antiprotons, and the difference between the two decreases with 
increasing energy. This effect can be understood either as a result of baryon stopping [5, 6] 
or by including the mean-field potentials in the hadronic phase, as shown in the AMPT 
model [7]. 

An intellectually more challenging experimental observation is the dependence of elliptic 
flow for 7r + and 7r~ as a function of the event-by-event fluctuating asymmetry parameter, 



where N + and A_ denote, respectively, the numbers of positively and negatively charged 
particles at a given acceptance region. According to the preliminary STAR data [8, 9], at 
a certain energy and centrality class, the elliptic flow for ir + is a decreasing function of A, 
and, for ir~ , the opposite trend is observed. In particular 



where (^2)^4 denotes an average elliptic flow at a given value of A, and c and r are parame- 
ters. The preliminary STAR data shows that in Au+Au mid-peripheral collisions at various 
energies r exp ~ 0.03 [8, 9], where both A and V2 are measured in \rj\ < 1. 

As shown in [10], this dependence cannot be understood in the baryon- (or isospin-) 
stopping scenario. Moreover, both UrQMD [11] and AMPT [12] models fail to describe the 
linear relation (1.2) [8]. Recently, an interesting interpretation of Eq. (1.2) was proposed 
as evidence of the effects of a chiral magnetic wave [13-15]. This effect is closely related to 
the chiral magnetic effect [16], discussed extensively in the recent literature, e.g., the recent 
review in [17]. 

In this paper, we demonstrate that Eq. (1.2) is consistent with the hydrodynamic picture 
of heavy- ion collisions with local charge conservation at freeze-out. We discuss the relation 
between Eq. (1.2) and the physics underlying the charge balance function [18], successfully 
described by hydrodynamic models, e.g., [19, 20]. In the following section, we derive a 
general relation between (v^)a and A, yielding a linear dependence. In section 3, we discuss 
two mechanisms that qualitatively and partly quantitatively reproduce the experimental 
data. In section 4, we detail state-of-the-art 3 + 1-dimensional hydrodynamic calculations 
with local charge conservation, observing a strong signal consistent with the experimental 
one, within a factor of 2. Section 5 gives our comments and conclusions. 

II. GENERAL RELATIONS 

In this section we derive the relation between elliptic flow at a given A, (v2) A , as a 
function of (1*2) ^ =0 an d A. Here, for simplicity, we assume that the distribution of A at a 




(1.2) 
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given centrality class is symmetric with respect to A = 0. 1 
Straightforwardly, we notice that for a sufficiently small A 



(N+ - N_) A = aA (N + + AL) 



A=0 



2a A (N+) 



A=0 ' 



(2-1) 



where a is a parameter whose value is close to unity 2 , see Eq. (1.1). Here, (.) A=0 denotes 
an average over all events characterized by A = 0. We readily notice that 



(N + + N-) A = 2 (N + ) A=0 (l - p 2 A 2 + ...) 



(2.2) 



The left-hand side of the above equation is invariant under the transformation A — > —A; 
thus, only the even powers of A are allowed. 3 These equations above support our evaluation 
of {N+) A and {N.) A 



(N ± ) A = (N + ) A=0 (l±aA- (3 2 A 2 + ...). 
To express elliptic flow at a given A as a function of A we write 

| / d 3 N + \ 

\ d(fd 2 p t I A=0 \ difd 2 p t ' 



I d 3 N+ \ / d 3 N^ 
\dipd 2 p t / A 



(2.3) 



(2.4) 



asym 



and analogously for negative particles. Here, (.) represents an average number of par- 
ticles responsible for a non-zero value of A. Integrating the above equation over (p and p t 
we obtain Eq. (2.3) that allows us to express (N + ) asym by (N + ) A=Q and A. Multiplying 
both sides of Eq. (2.4) by cos(2(^) and integrating over an available phase space, assuming 
A << 1 such that (1 ± A)' 1 m 1 =p A and neglecting the A 2 terms, we obtain 



A=0 



- aA 



A=0 



and 



A=0 



+ aA 



A=0 



- {V 



asym 



asym 



(2.5) 



(2.6) 



where {v 2 )a denotes elliptic flow at a given value of A, and {v 2 ) aS ym represents an average 
elliptic flow of particles responsible for a non-zero value of A. The elliptic flow for all events 
with A = is denoted by (u^) A=0 . 

As is evident from the above equations, a linear dependence on A is expected naturally 
only if the elliptic flow of particles contributing to a non-zero asymmetry differs from the 
elliptic flow of particles leading to A = 0. In the next section, we discuss two mechanisms 
that naturally lead to (v2)a=o > (1^2) asym ; being qualitatively consistent with the preliminary 
data. 



1 So that we can expand all variables around A = 0. 

2 If the same particles are used to calculate A and {N±) A , then a ~ 1. In the STAR measurement, (N±) A 
and A are calculated from similar sets of particles, but not identical ones. This distinction may engender 
a value that differs slightly from 1. 

3 To obtain A = 0.1, for example, we could increase N+ — iV_, or/and decrease N+ + N-. Thus, for larger 
values of A, N + + N- is expected to decline as a function of A. Wc verified this in our simple model 
discussed in the next section, where f3 fa 2.3, and a = 0.96. 
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III. TWO MECHANISMS 



In this section, we discuss two mechanisms that offer qualitative and partly quantitative 
agreement with the preliminary STAR data. 

From various measurements of the balance function in rapidity, we know that electric 
charges balance each other with a characteristic distance between the charges of about 0.5 
unit of rapidity [21]. This finding is confirmed successfully in the hydrodynamic model with 
locally conserved charge production at freeze-out [19, 20]. 4 Consequently, a pair of pions 
with opposite charges originating from a fluid element or a resonance, create a non-zero A 
in a given rapidity window provided that one pion is produced inside of this window and 
the other outside it. If both charges are produced inside a bin they do not contribute to 
N + — AL. In the STAR experiment, the elliptic flow is measured in \r)\ < 1 that is broader 
than the typical width of the balance function. 

In Fig. 1, we schematically depict the first mechanism, leading to Eq. (1.2). Fluid 
elements and resonances, in short clusters, located close to rj = usually produce both 
pions inside the pseudorapidity window. Only clusters located close to \rj\ = 1 typically 
create a non-zero asymmetry parameter. From the PHOBOS data [22], we know that 
elliptic flow of charged particles as a function of pseudorapidity can be roughly described as 
V2{v) ~ 1 — M/6 (further discussed in section 5). Thus, clusters located close to |?7| = 1 have 
smaller elliptic flow compared to clusters located close to r\ — 0. Consequently, the elliptic 
flow of pairs responsible for a non-zero asymmetry is smaller than that of pairs leading to 
A = 0. Using Eq. (1.2) we conclude that {v%)a decreases and {v^)a increases as a function 
of A. 5 



-1 6 i i 

FIG. 1: Clusters (resonances, fluid elements), denoted by the blue dots, positioned close to r/ = 0, 
decay into final particles located inside the rapidity bin \rj\ < 1 and clearly cannot create a non-zero 
asymmetry, A. Only clusters located close to the boundary can produce one pion inside and one 
pion outside of our rapidity bin. From the PHOBOS data, we know that elliptic flow is slightly 
smaller away from r] = 0. Consequently, the elliptic flow of particles creating the asymmetry, 
A ^ 0, is smaller than that of particles leading to zero asymmetry, A = 0. 

The second mechanism is explicitly related to the late locally conserved charge produc- 
tion in the hydrodynamic picture, see Fig. 2. Clusters with small transverse momentum 



4 In contrast to models with an initial charge creation (such as UrQMD) that cannot describe the experi- 
mental data on the balance function [21]. 

5 By selecting events with A > 0, for example, we trigger configurations wherein positive pions are inside 
a rapidity window and the negative ones are outside it. In this way, we add positive particles with small 
i>2 and, effectively, the i>2 for positive particles is lowered. Since the number of negative particles with a 
small vi is reduced, thus effectively, their vi is increased. 
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(resonances) or small transverse velocity (fluid elements) usually decay into pairs of pions 
with a large rapidity separation (a broad balance function). In contrast, clusters with large 
transverse momentum or velocity engender a small rapidity separation (a narrow balance 
function). Undoubtedly, clusters with a broad balance function are more likely to produce 
a non-zero charge asymmetry parameter, A. On the other hand, we know that the elliptic 
flow of clusters with a small transverse momentum (velocity) is less than that of clusters 
with a high transverse momentum, e.g., [3]. Consequently (^2) asym is smaller than (v 2 ) A=0 , 
so leading to the linear dependence of (v^)a — ( v £)a on A. 

- + 
- + 

3 o i "n 

FIG. 2: Clusters (resonances, fluid elements), denoted by the blue dots, with low transverse mo- 
mentum (velocity) produce particles separated more in rapidity compared to clusters with high 
transverse momentum. Consequently, non-zero values of A are mainly driven by clusters with a 
small pt that are characterized by a small vi. Thus, the elliptic flow of particles creating asymmetry, 
A ^ 0, is smaller than that of particles leading to zero asymmetry, A = 0. 

To quantify these effects, we undertook a simplified Monte Carlo calculation. 6 As shown 
in Fig. 3, the signal we obtained is of the same order-of-magnitude as the data. In fact, we 
obtained r « 0.02, see Eq. (1.2), that is only a factor of 1.5 smaller than the preliminary 
STAR data, r exp « 0.03. In our calculation, both A and v 2 are calculated for particles in 
\r}\ < 1. However, we note that our results represent only an order-of-magnitude estimate, 
and more advanced calculations are required to draw definite conclusions. 

We checked that the effect schematically presented in Fig. 1 is responsible approximately 
for 35% of the signal shown in Fig. 3. We also checked that calculating v 2 in \rj\ < 0.2 and A 
in \t}\ < 1 reduces our final results by about 10%. However, when calculating v 2 in \rj\ < 1.5, 
and A in \rj\ < 1, we obtained a significantly smaller value, r ps 0.006. Taking v 2 in \rj\ < 2 
and A in \r]\ < 1 lowers our results approximately by a factor of 10. Finally, we mention 



6 In each event, we sampled clusters from the Poisson distribution with an average number of clusters, 
(N c ) = 570. All clusters were uniformly distributed in pscudorapidity over < 3 and in transverse 
momentum according to the thermal distribution, e~ 2p *^ Pt ' > with (p t ) = 0.7 GeV. A cluster with a given r\ 
and pt decays into two oppositely charged particles, located symmetrically around r\. The pscudorapidity 
distance between them is sampled from the Gaussian distribution with its standard deviation a = 1/(1 + 
Pt / {Pt}), so that for a cluster with p t = (p t ), a — 0.5. This dependence on p t approximately agrees with 
a rho meson decay at a given p t into two pions. Finally, we sampled the azimuthal angles of the two 
particles using 1 + 2v2(pt,v) cos(2y)) with v 2 (pt,r?) given by 0.04(1 — |?7| /Q)pt/ (pt) for p t < 2 GeV, and 
0.04(1 — \r)\ /6)2/ (pt) for p t > 2 GeV. A and i>2 are calculated for particles in \rj\ < 1. In our calculation 
we obtained (V2) « 0.036, and approximately 380 charged particles in the midrapidity region \r]\ < 1, 
approximately agreeing with 200 GeV Au+Au collisions in the 30 — 40% centrality class. 
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FIG. 3: Our results for the elliptic flow coefficient for positive and negative particles as a function 
of asymmetry parameter A = (N + — N_)/(N + + iV_). In our simplified model we obtain r ~ 0.02, 
see Eq. (1.2), to be compared with the preliminary STAR data r cxp ~ 0.03. Note that our results 
represent only an order-of-magnitude estimate. 



that in our simplified analysis we neglected experimental cuts in the transverse momentum 
which may slightly modify our numbers. 

It would be desirable to perform a more realistic calculations in an event-by-event 3+1- 
dimensional hydrodynamic model, wherein the effect presented in Fig. 2 is naturally present. 
This problem is discussed in the next section. 



IV. HYDRODYNAMIC CALCULATION 



The collective flow in heavy-ion collisions can be reproduced satisfactory within relativis- 
tic hydrodynamics [1-3]. Statistical emission and resonance decay at freeze-out account 
for part of the observed charge balancing correlations. Additionally, local charge balancing 
in particle production [18] can be included in the event generator [20], yielding the right 
qualitative description of the one- and two-dimensional charge dependent correlations. We 
used a 3 + 1-dimensional hydrodynamic model with statistical emission at freeze-out includ- 
ing the local charge conservation mechanism and resonance decays [23]. In order to study 
non-flow correlations, particle emission from the freeze-out hypersurface should be done via 
Monte-Carlo event generation. The default implementation of the statistical emission in 
THERMINATOR [24] generates primordial particles with the Poisson distribution for each 
type of particle. In the present calculation, we fix the total charge produced in the fire- 
ball, for simplicity fixing the multiplicities for all particles emitted from the fireball in each 
event. The charge asymmetry measured in a limited acceptance window mainly is due to 
the separation of an unlike charged pair as described in the previous section. 

The statistical emission code [24] follows resonance decays, and such charge correlations 
are realistically described in the model. Together with the effects of local charge conserva- 
tion, correlating the emission of primordial particles, it yields a good quantitative description 
of charge balance functions in semi-central Au-Au collisions [23]. To study the weak charge 
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FIG. 4: The charge asymmetry dependence of tt + and ir~ elliptic flow coefficients in the hy- 
drodynamic model followed by statistical emission with local charge conservation. We obtained 
r = 0.012 i 0.004 compared with the preliminary STAR data, r cxp « 0.03. 



splitting signal, a large statistics is required. We generated 10 6 events for the centrality 
30-40% wherein the measured signal is the strongest, and where the hydrodynamic model 
agrees best with the data on charge balance functions. The elliptic flow, V2, for 7r + and 7r~ 
is calculated for events with different charge asymmetry in the interval \r]\ < 1, using two 
subevents with a pseudorapidity separation | Ar)\ > 0.3 to define V2, following the experimen- 
tal procedure [8]. The Monte-Carlo calculation with local charge conservation and resonance 
decays exhibits the expected signal, viz., a decrease of V2(tt + ) and an increase of V2(tt~) with 
the charge asymmetry A in the event (Fig. 4). The slope parameter, r ~ 0.012, is about one 
half of the observed magnitude. We note that the pseudorapidity dependence of the elliptic 
flow ^2(77) observed experimentally [22] is noticeably underestimated in 3 + 1-dimensional 
hydrodynamic simulations [23, 25]. The steeper pseudorapidity dependence of v 2 is expected 
to improve the agreement between the calculated charge splitting of the elliptic flow with 
the experimental findings. In the next section, we comment further on this issue. 

Finally, we discuss the centrality dependence of the signal. Two effects that are not 
well controlled in hydrodynamic models may cause the decrease of the slope parameter for 
peripheral events. First, the balance function in pseudorapidity is widening for peripheral 
events. Second, for peripheral events, correlations due to elliptic flow are diluted by the larger 
contribution of particles emitted from the nonthermal corona in the interaction region, while 
the charge asymmetry acquires contributions both from the core and the corona. 



V. COMMENTS AND CONCLUSIONS 



Several comments are warranted. 

(i) In this paper we discussed two mechanisms leading to elliptic flow splitting as a 
function of the asymmetry parameter, A. We note that we expect a similar effect for higher 
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harmonics, v n . In particular, for the triangular flow 




although the signal is expected to be significantly weaker in comparison to the case for v 2 - 
We expect that r, see Eq. (1.2), should be smaller roughly by a factor of (v 2 ) / (V3) ~ 3 [26]. 
This problem will be discussed elsewhere. 

(ii) The key ingredient in our analysis is the assumption of local charge conservation 
at freeze-out. In this scenario, particle pairs with higher momentum are more strongly 
collimated in rapidity than are pairs with smaller momentum, see Fig. 2. This allows us to 
understand the centrality dependence of the rapidity balance function [21], in contrast to 
models with initial charge creation (such as UrQMD). Consequently, in models without late 
local charge conservation we do not expect to observe the elliptic flow splitting discussed in 
this paper. 

(iii) According to the PHOBOS data, the elliptic flow of charged particles in the midrapid- 
ity region changes significantly as a function of pseudorapidity [22]. However, the results of 
the STAR Collaboration indicate that elliptic flow in the midrapidity region weakly depends 
on r\ [27], so being consistent with hydrodynamic calculations (see section 4). Consequently, 
the mechanism presented in Fig. 1 may be suppressed if we take the STAR data into ac- 
count. However, this is not a problem because, as we discussed in section 3, this mechanism 
is responsible for only 25% of the measured signal, and 35% of the signal presented in Fig. 3. 

In conclusion, we studied the dependence of the elliptic flow coefficients for positive and 
negative particles as a function of the event-by-event charge asymmetry parameter. Recently, 
this phenomenon was interpreted as evidence of the chiral magnetic wave. In this paper, we 
argued that the origin of this effect may be less exotic and actually is consistent with the 
hydrodynamic picture of heavy-ion collisions with late local charge conservation. 

We argued that particle pairs leading to a non-zero asymmetry parameter, A ^ 0, are 
characterized by a smaller elliptic flow in contrast to particle pairs resulting in A = 0. We 
showed that this fact alone is sufficient to qualitatively explain the preliminary STAR data, 
in particular, a linear dependence of (v 2 )a — (v 2 )a, as a function of A. Our quantitative re- 
sults, based on a simplified Monte Carlo estimation and the state of the art 3 + 1-dimensional 
hydrodynamic model calculations, are in agreement, within a factor of 2, with the prelimi- 
nary STAR data. 
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